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Abstract 

Single neutral pion production via muon antineutrino charged-current interactions in plastic scintillator (CH) is studied 
using the MINERvA detector exposed to the NuMI low-energy, wideband antineutrino beam at Fermilab. Measurement 
of this process constrains models of neutral pion production in nuclei, which is important because the neutral-current 
analog is a background for P e appearance oscillation experiments. The differential cross sections for 7 r° momentum and 
production angle, for events with a single observed 7 r° and no charged pions, are presented and compared to model 
predictions. These results comprise the first measurement of the 7 r° kinematics for this process. 
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1. Introduction 

Neutrino- and antineutrino-induced interactions at en¬ 
ergies of a few GeV are a proving ground for weak in¬ 
teraction phenomenology in nuclei [1]. Measurements in 
this energy range are important because neutrino oscil¬ 
lation experiments [2, 3] need detailed understanding of 
the large variety of processes allowed. As a result, there 
is a growing body of new high-quality measurements for 
neutrino-nucleus interactions [4-7]. In particular, neutrino 
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charged-current neutral pion production has become an 
important benchmark providing new challenges to theo¬ 
ries describing this process [8-11]. 

Most of the published data for pion production in nuclei 
uses neutrino beams. Neutral pion production in nuclei for 
anti-neutrinos, however, is much less studied. Only one 
data point for this channel in the few-GeV energy range 
exists in the literature: a measurement of u^p —> p, + mr° 
from SKAT in a heavy liquid (CF 3 Br) bubble chamber 
based on 20 events at an average neutrino energy of 7 
GeV [12]. Measurements of tt° production by neutrinos 
have been made in deuterium bubble chambers [13-16] 
and more recently on nuclear targets in the 0.1-1 GeV 
energy range using the MiniBooNE detector [17] with a 
mineral oil (CH 2 ) target, and using the SciBar detector in 
K2K and SciBooNE experiments [18, 19] with plastic scin¬ 
tillator (CH). These recent measurements, as well as data 
on charged pion production [7, 20] and neutral pion pro¬ 
duction [21] have been difficult for event generators and 
theoretical calculations to describe accurately [8-11]. Ev¬ 
ery prediction must have a model for tt° production from 
nucleons; all use isospin decompositions within a helicity 
formalism that are tuned to available data [22]. 

Charged-current single 7r° (l7r°) production in the few- 
GeV region is modeled both as decays of nucleon reso¬ 
nances (most strongly the A(1232)) and nonresonant pro¬ 
cesses. The production in nuclei can be either direct, 
through the reaction v^p —> or indirect, for ex¬ 

ample, through charge exchange (CEX) of a charged pion 
in the nucleus, p~K~ —> mr° or n7r + —> p tt°. New data will 
provide useful tests of both neutrino-induced resonance 
production and final state interaction (FSI) models. 

Neutrino interaction measurements are also important 
to the analysis of neutrino oscillation experiments [2, 3, 
23, 24]. These experiments require the neutrino flavor 
be identified and the neutrino energy to be reconstructed 
on an event-by-event basis. An accurate modeling of 
these particles requires knowledge of both the underlying 
neutrino-nucleon interactions and of the final-state modifi¬ 
cations that arise within the target nuclei (such as carbon) 
of which the massive oscillation detectors are comprised. 
Pion production is a source of backgrounds and systematic 
uncertainties in neutrino oscillation experiments. Neutral- 
current 7T° production, for example, is a dominant back¬ 
ground in v e ( v e ) appearance experiments because the 7r° 
can mimic a final state electron (positron). In addition, 
experiments that reconstruct neutrino energy by identify¬ 
ing quasielastic events, veN(n) —> £~p or vgN{p) —> £ + n, 
interactions in which a pion is produced but then absorbed 
in the target nucleus can be mistaken for quasielastic sig¬ 
nal and yield an incorrect estimate of the incident neutrino 
energy. 

New measurements of l7r° production by charged- 
current interactions in plastic scintillator (CH) using 
the MINERvA detector are presented. Flux-integrated 
single differential cross sections as a function of 7r° mo¬ 
mentum and production angle for events with a single ob¬ 


served 7r° and no exiting the interaction nucleus have 
been measured and are compared to predictions from the 
GENIE [25], NuWro [26, 27], and NEUT [28] event gener¬ 
ators. 

2. Experiment 

The data presented here were taken using the MINERvA 
detector and the wideband antineutrino beam produced 
by the NuMI beamline in the low-energy mode [29] with 
a mean energy of 3.6 GeV. The antineutrino flux is esti¬ 
mated from a simulation of the neutrino beamline based 
on Geant4 [30, 31], with hadron production in the simula¬ 
tion constrained by proton-carbon external data [32-34]. 
The MINERvA detector consists of a fully active region of 
scintillator strips surrounded on the sides and downstream 
end by electromagnetic and hadronic calorimeters, and is 
described in detail in Ref. [35]. There are three orienta¬ 
tions (views) for the strips (X,U,V), offset by 60° from 
each other, which enable three-dimensional reconstruction 
of particle trajectories. The X view is sampled twice as 
often as the other views. The downstream edge of the 
MINERvA detector is located 2 m upstream of the MI¬ 
NOS Near Detector, a magnetized iron spectrometer [36] 
used in this analysis to reconstruct the momentum and 
charge of muons. The transport of particles from neutrino 
interactions in the detector is simulated by a Geant4-based 
program. The readout simulation is tuned so that both the 
photostatistics and the reconstructed energy deposited by 
momentum-analyzed through-going muons agree between 
data and the simulation. The detector simulation of single 
particle responses is validated using testbeam data taken 
with a scaled-down version of the MINERvA detector [37]. 

Neutrino interactions are simulated using the GE¬ 
NIE 2.6.2 neutrino event generator. Details concerning 
GENIE and its associated parameters are described in 
Ref. [25]. For baryon resonance production, the formal¬ 
ism of Rein-Sehgal [22] is used with modern resonance 
properties [38]. Non-resonant pion production is simu¬ 
lated using the Bodek-Yang model [39] and is constrained 
below IV = 1.7 GeV by neutrino-deuterium bubble cham¬ 
ber data [40, 41]. Pion and nucleon FSI are modeled in 
GENIE using a parameterized intranuclear cascade model, 
with the full cascade being represented by a single inter¬ 
action. For all models fitting data for light nuclei such as 
carbon, the pion most often has one interaction as it prop¬ 
agates through the nucleus. For each interaction, choice 
of the channel (e.g. charge exchange) is based on total 
cross section data and calculations [42, 43] and the kine¬ 
matics and multiplicity of the final state are taken from fits 
to more detailed data. At the energies important for this 
measurement, the hadron-nucleus cross sections have large 
uncertainties. For example, the cross section of 7r“ —> 7r° 
CEX on carbon has an uncertainty of about 50%. Re¬ 
action cross sections for 7r° are estimated from measured 
cross sections for scattering using isospin symmetry 
or theoretical models. The most significant advantages of 
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the single interaction used in GENIE’s FSI model are the 
ability to exactly reweight and to characterize each event 
with a single final state channel. The uncertainties in the 
FSI model are evaluated by varying its strength within 
previously measured hadron-nucleus cross-section uncer¬ 
tainties. 

This analysis uses data taken between October 2009 and 
February 2012 with 2.01 x 10 20 protons on target (POT) 
in the mode. About half of the exposure (0.945 x 10 20 
POT) was taken during construction with the downstream 
half of the detector. In this period the ArgoNeuT detec¬ 
tor [44] was situated between the MINERvA and MINOS 
Near detectors. Because the two sub-samples of the data 
have different efficiencies, they were analyzed separately 
and their results combined. 

3. Event reconstruction and selection 

The MINERvA detector records the charge and time of 
energy depositions (hits) in each scintillator strip. Hits 
are first grouped in time and then clusters of energy are 
formed by spatially grouping adjacent hits in each scintilla¬ 
tor plane. Clusters with energy more than 1 MeV are then 
matched among the three views to create a track. The per- 
plane position resolution is 2.7 mm and the angular resolu¬ 
tion of the muon track is better than 10 mrad [35] in each 
view. The /i + is identified by matching a track that ex¬ 
its the back of MINERvA with a positively-charged track 
entering the front of MINOS. The reconstruction of the 
muon in the MINOS spectrometer gives a typical momen¬ 
tum resolution of 11%. Event pile-up causes a decrease 
in the muon track reconstruction efficiency. This effect 
is studied in both MINERvA and MINOS by projecting 
tracks found in one of the detectors to the other and mea¬ 
suring the mis-reconstruction rate. This results in a —4.4% 
(—1.1%) correction to the simulated efficiency for muons 
below (above) 3 GeV/c. 

The event vertex, defined as the most upstream cluster 
on the muon track, is restricted to be within the central 
108 planes of the scintillator tracking region and no closer 
than 22 cm to any edge of the planes. These requirements 
define a fiducial volume with a mass of 5.47 metric tons. 
Due to the requirement that the /i + is tracked in MINOS 
for charge and momentum measurement, the detection ef¬ 
ficiency has a strong dependence on muon angle and 
momentum \p^\, which drops to zero for events with 6 M 
greater than 20 degrees or \p tl | less than 1.0 GeV/c. The 
corrections for the angle and momentum efficiency are es¬ 
timated from the event simulation. Only events with a 
single track at the vertex, the /i + matched to MINOS, 
are used in order to reject events including charged pion 
production. Accepted events are passed to the tt° recon¬ 
struction. 

The neutral pion has a lifetime of 8.52 x 10~ 17 s and de¬ 
cays into two photons with a branching ratio of 98.8% [38], 
so the two photons appear to come from the event vertex. 
Plastic scintillator has a radiation length X 0 ~ 40 cm, 


which allows the two photons to convert by e + e~ pair 
production or Compton scattering far away from the ver¬ 
tex, thus producing isolated energy deposits. The photons 
can also have significant energy leakage or escape the de¬ 
tector without conversion. Furthermore, energy deposits 
produced by neutrons from the same neutrino interaction 
can be mistaken for low-energy photon conversions, fur¬ 
ther complicating the pattern of energy deposits. The 7 r° 
reconstruction must correctly group these energy deposits 
into the two photons. 

The reconstruction can be separated into two steps: pat¬ 
tern recognition and kinematic reconstruction. The pat¬ 
tern recognition builds upon the knowledge of the vertex 
location of the event. In the X view, clusters that are close 
in polar angle with respect to the vertex but can be sepa¬ 
rated in radial distance from the vertex are grouped into 
photon candidates. Then, for each candidate, clusters in 
the U and V views consistent with the stereo condition are 
added. Photon candidates must have clusters from at least 
two views for three-dimensional direction reconstruction. 

In the second step, photon position, direction, and en¬ 
ergy are determined from the clusters that have been as¬ 
signed to the photons. The photon direction is recon¬ 
structed from the cluster energy-weighted slopes in each 
view. The photon vertex is defined by the closest cluster to 
the event vertex on the photon direction axis. The photon 
energies are reconstructed by calorimetry using calibration 
constants determined from the simulation. The overall cal¬ 
ibration constant that sets the absolute energy scale is de¬ 
termined by matching the peak in the 77 invariant mass 
distribution to the n° nominal mass of 134.97 MeV/c 2 [38]. 
This procedure is done separately for data and simulation 
which enables correction for a difference in energy scales 
of 5% between the data and simulation. Finally, the 7 r° 
momentum is calculated from momentum conservation, 
p n 0 = ki+k 2 , where fe; are reconstructed photon momenta. 
The 7 T° is reconstructed with a 25% energy resolution and 
3.5 degrees angular resolution in each view. The incom¬ 
ing neutrino energy is reconstructed from the p + and 7 r° 
4-moment a using 

E r v ec = E^ + E^o+Tn (1) 

rr _ 1 K-£'A~" f vi) + ( E K° — p![o)] 2 + {p^ +P^o) 2 

n II || 1 

z m N - (E^ - pji) - ( E n o - p'l 0 ) 

where are the transverse and longitudinal compo¬ 

nents of momentum, respectively. It is assumed that the 
initial nucleon is at rest and that the 7 r° is produced to¬ 
gether with a nucleon. The neutrino energy is recon¬ 
structed with 10% resolution. The 77 invariant mass m 77 
is reconstructed from the photon energies E\^ and the 
separation angle 9 11 between the two photons using 

m 2 7 = 2E 1 E 2 {1 - cos0 77 ). ( 2 ) 

The two reconstructed photons are required to convert 
at least 15 cm (0.36X 0 ) away from the vertex to further re¬ 
ject charged-pion backgrounds. In addition, it is required 
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that E^ ec is between 1.5 and 20 GeV. The lower energy cut 
is needed due to MINOS acceptance while the upper cut 
is to reduce flux uncertainties. Finally, it is required that 
m 77 is between 75MeV/c 2 and 195MeV/c 2 . The selected 
sample has 1304 events. The total selection efficiency is 6% 
and purity 55%, according to the simulation. The back¬ 
ground is dominated (70% of the total background) by 
events with at least one 7r° in the detector. Half of this is 
due to multi-pion production, 7T°+7r , where the is not 
tracked, while the other half has a secondary 7r° produced 
by 7r~ —» 7r° CEX or nucleon scattering in the detector 
but outside the primary interaction nucleus. The non-7T° 
background is mostly due to energy deposits by 7r - and 
neutrons which are mistakenly identified as photons, and 
accounts for the remaining 30% of the total background. 

Figure 1 shows the m 77 distributions for both data and 
simulation of the selected sample before the invariant mass 
cut. There is a clear peak centered around the n° nom¬ 
inal mass in both data and simulation. The distribution 
from simulation is broken down into signal and background 
components. The signal is defined as antineutrino charged- 
current events with single 7r° and no 77^ escaping the nu¬ 
cleus. The background is anything else that is not signal. 
By this definition, it is possible for signal events to have the 
7r°s mis-reconstructed from non-7r° energy deposits. The 
signal events at high ?n 77 , outside the signal mass width, 
have one or both candidate photons reconstructed from 
neutron energy deposits. The same mis-reconstruction 
happens to signal events in the selected sample but at a 
smaller rate (14%). The 7r° momentum and angular shapes 
of these events are found to be similar to the rest of the 
signal events. The enhancement in the background distri¬ 
bution around the 7 r° mass is due to background events 
with at least one 7r° in the detector. 

After event selection, the selected sample still has sub¬ 
stantial background to be subtracted statistically from 
the total distribution for each observable. The back¬ 
ground distribution for each observable is estimated from 
the simulation with the total background rate constrained 
by data. This significantly reduces the uncertainties in 
the estimated background. The background rate is de¬ 
termined from a binned extended maximum likelihood fit 
of an invariant mass model to the data [45]. The invari¬ 
ant mass model M(m 77 ) is a two-component model con¬ 
structed from the m 77 distributions of signal and back¬ 
ground events, 

M(m 77 ) N s ig^^-sig (t7T 77 ) T Nbkg^^bkg(,TH"y r y') : (3) 

where M si9 (m 77 ), Mbkgim-yj) are the shapes of the sig¬ 
nal and background m 77 distributions from the simula¬ 
tion, respectively. The expected numbers of signal and 
background events N sig ,Nbk g in the range 0-500 MeV/c 2 
are the parameters determined from the fit. After the fit, 
a background rate of 541±37 events is obtained by in¬ 
tegrating the curve Nbk g Mbk g {m 11 ) over the mass peak 
region from 75MeV/c 2 to 195MeV/c 2 , the same range as 


required by the event selection. The fit reduces the back¬ 
ground normalization by a factor of 0.8 compared to the 
simulation prediction. 



m yy (MeV/c 2 ) 


Figure 1: Distribution of the invariant mass of the 77 pair. Data 
are shown as solid circles with statistical error bars. The shaded his¬ 
tograms show the Monte Carlo predictions for CC1 7T° signal (top) 
and for background (bottom). The signal is defined as antineutrino 
charged-current events with single 7r° and no 77 + escaping the nu¬ 
cleus. The background is anything else that is not signal. The ver¬ 
tical lines indicate the invariant mass cut, 75MeV/c 2 < m 77 < 
195 MeV/c 2 . 

The estimated background is subtracted from the to¬ 
tal distribution. This background-subtracted data distri¬ 
bution is unfolded to account for detector resolution ef¬ 
fects using a Bayesian unfolding method [46] with two 
iterations. The migration matrix used in the unfolding 
is obtained from the simulation. The reconstruction ef¬ 
ficiency and acceptance corrections are made by dividing 
by an efficiency curve derived from the simulation. Di¬ 
viding this corrected distribution by the integrated flux 
(2.5 x 10 _8 il /i /cm 2 /POT) of antineutrinos with energies 
in the range 1.5-20 GeV, and the number of nucleons in 
the fiducial volume (3.30 x 10 30 ) gives the differential cross 
section. 

The total uncertainties on the measured cross sections 
are 20-25% with comparable contributions from statisti¬ 
cal and systematic uncertainties. The largest systematic 
uncertainty is due to the integrated flux uncertainty at 
10%. A detailed discussion of the flux systematic uncer¬ 
tainty is presented in Ref. [6]. The next largest contribu¬ 
tion to the systematic uncertainty is the 8% normalization 
uncertainty from the background fit. The other smaller 
contributions include the neutrino cross-section models, 
FSI models, and detector simulation systematics. These 
systematic uncertainties enter the measured cross sections 
through background subtraction, detector resolution, and 
efficiency corrections. The neutrino cross-section model 
and FSI model uncertainties are evaluated by GENIE. The 
systematic uncertainties on the estimated background is 
small since the background rate is constrained by data. 
Systematic errors in the detector response must be inde- 
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pendently estimated. The uncertainty in neutron response 
is evaluated by changing the neutron inelastic cross section 
within experimental uncertainties [47-54] through event 
reweighting. The reweighting is applied to the leading 
neutron in the event. A large fraction of the secondary 
7T° in the background is estimated to arise from 7r“ —> 7 r° 
CEX, for which the cross sections are poorly known. The 
effect of this uncertainty on our measurement is evalu¬ 
ated by changing the CEX cross section within its uncer¬ 
tainty of ±50% [43, 55, 56], and then re-measuring the 
cross sections. Finally, the electromagnetic energy scale 
uncertainty (2.2%) is taken from the fitted mean uncer¬ 
tainty of the data m T7 distribution. Flux-integrated single 
differential cross sections in x° momentum and angle with 
statistical, systematic, and total uncertainties are summa¬ 
rized in Tables. 5 


4. Results and discussion 

The measured differential cross sections as function of x° 
momentum and angle with respect to the beam direction 
are shown in Figures. 2 and 3, respectively. The data are 
compared to the predictions from GENIE with and with¬ 
out FSI. Above 0.7 GeV/c, FSI effects have little influence 
on the 7T° momentum distribution, and both predictions 
are in good agreement with the data. For momenta below 
0.3 GeV/c, inclusion of FSI gives an increased and shifted 
cross section relative to the no FSI case. This trend is ex¬ 
hibited by the data; GENIE calculations with and without 
FSI give a x 2 of 25.4 and 50.0 for 11 degrees of freedom 
(dof), respectively. For the distribution of 7r° production 
angle in Fig. 3, inclusion of FSI into the GENIE simula¬ 
tion results in a mild flattening of the distribution with 
no significant improvement in the x 2 compared to the no 
FSI case, 16.7 versus 16.0 for 11 degrees of freedom. Thus 
the effects of FSI are more pronounced with respect to 
7T° momenta than with 7r° production angle. This situa¬ 
tion likely reflects the influence of the A(1232) resonance, 
which gives a particularly strong momentum dependence 
to the pion-nucleus interaction for p * ~ 0.26 GeV/c where 
the pion-carbon cross section is maximum. 

Figures. 2 and 3 also show predictions including FSI 
from the NuWro and NEUT event generators. Any pre¬ 
diction requires knowledge of v^N —> p + x°N reactions. 
Because of the dearth of data for these channels, the calcu¬ 
lations use isospin relations to extrapolate from other pion 
production channels [22, 57]. Both NuWro and NEUT 
use full cascade models [58] for the FSI description. The 
X 2 /dof values for the NuWro and NEUT comparisons are 
25.0/11 and 24.9/11 for the 7r° momentum and 11.8/11 
and 14.0/11 for the 7r° production angle, respectively. 
These \ 2 values indicate that a common level of agreement 
is achieved by all three generators for x° momentum, and 
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Figure 2: Differential cross section for l7T° production as function of 
7 r C) momentum. Data are shown as solid circles. The inner (outer) 
error bars correspond to statistical (total) uncertainties. The solid 
(dashed) histograms are GENIE prediction with (without) FSI, the 
long-dashed histogram is the prediction from the NuWro generator, 
and the dot-dashed histogram is the prediction from NEUT. 



6*0 (deg.) 


Figure 3: Differential cross section for I 7T° production as function of 
the 7r° polar angle. Data are shown as solid circles. The inner (outer) 
error bars correspond to statistical (total) uncertainties. The solid 
(dashed) histograms are GENIE prediction with (without) FSI, the 
long-dashed histogram is the prediction from the NuWro generator, 
and the dot-dashed histogram is the prediction from NEUT. 


that modest differences with predictions versus data are 
observed for 7r° production angle. 

There are uncertainties in the FSI used in all calcula¬ 
tions. Each assumes the pion FSI after production in the 
nuclear medium is the same as for pion beams; this as¬ 
sumes no medium effects beyond Fermi momentum and a 
simple binding energy. Since x° beams are not possible, 
isospin relations must be used for x° FSI. Experiments 
like these are valuable for testing these approximations. 
In spite of these uncertainties, the calculations give ade¬ 
quate descriptions of these new data. 
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Figure 4: Same cross section data as Fig. 2. The stacked histograms 
show a decomposition of the l7r° signal into different FSI channels as 
predicted by GENIE. Description of the FSI channels follows (from 
top to bottom): 1) Multi-7r —>■ 7r°: multi-pion produced by the pri¬ 
mary interaction and all other pions are re-absorbed inside the nu¬ 
cleus, except a 7T°, 2) tv~ —>■ 7r°: a 7t _ produced by the primary in¬ 
teraction, then charge exchanges inside the nucleus, 3) 7r° produced 
by the primary interaction and then undergoing inelastic scattering 
inside the nucleus, 4) 7r° produced by the primary interaction and 
then undergoing elastic scattering inside the nucleus, and 5) 7r° pro¬ 
duced by the primary interaction and exiting the nucleus without 
interacting. 

Figures. 4 and 5 show decompositions of the 7r° momen¬ 
tum and angle spectra of Figures. 2 and 3 according to the 
FSI channels as predicted by the GENIE simulation, al¬ 
lowing for a more detailed interpretation. This shows how 
each FSI process changes the spectrum and gives an indi¬ 
cation of what refinements are needed for better agreement 
with data. GENIE predicts that about 80% of the events 
undergo some FSI. In the momentum spectrum, absorp¬ 
tion events (not shown in the spectrum because the pion 
disappears in the nucleus) preferentially deplete the region 
centered at p n « 0.26 GeV/c, the momentum where the 
pion-carbon total reaction cross section is a maximum [42]. 
On the other hand, pion inelastic scattering and CEX 
interactions shift the pion momentum from high values 
to lower values, therefore contributing to the buildup of 
events at p n « 0.2 GeV/c seen in the data. The contribu¬ 
tions from multi-pion production followed by absorption, 
elastic scattering, and the non-interacting component have 
no noticeable effect on shape. 

For the angle spectrum FSI decomposition (Fig. 5) the 
effects are very different. The interactions displayed in 
the figure have a small influence on shape and the angular 
distribution of the A(1232) decay becomes important. The 
inelastic and CEX interactions are largely responsible for 
the buildup of events at backward angles. 

Comparison of the observed spectral shape in Fig. 4 to 
the GENIE prediction suggests that an increase in inelastic 
scattering together with compensating reduction in elastic 
and/or multiple-pion production with absorption, could 


Figure 5: Same cross section data as Fig. 3 and same decomposition 
as with Fig. 4. The stacked histograms show a decomposition of the 
l7r° signal into different FSI channels as predicted by GENIE. 

improve the agreement with data. On the other hand, 
Fig. 5 indicates such changes would worsen the agreement 
with the data for pion production at backwards angles. 
Thus it appears that a refined description might require 
separate, independent adjustments to the two spectra. 

While the previously reported information on the reac¬ 
tion studied here is too limited to enable comparisons, it 
is useful to compare with the recently reported MINERvA 
observations on v M induced charged pion production [20]. 
The latter data was obtained using a low-energy exposure 
in the same beamline, and the cross section normaliza¬ 
tions are carried out in a similar way for both analyses. 
Existence of both results provides stronger constraints on 
the calculations. The 7r° momentum range in this analy¬ 
sis is wider than the range shown for the charged pions of 
Ref. [20] (0-1.5 GeV/c versus 0.1-0.5 GeV/c) because the 
7r + data is limited by tracking and particle identification 
requirements. Nevertheless, both analyses show that GE¬ 
NIE predictions are significantly improved when FSI are 
accounted for. Both results show a peak in the momen¬ 
tum distribution at p n ~ 0.2 GeV/c which is seen in the 
GENIE calculations, but do not have the correct distribu¬ 
tion at energies near this peak. The GENIE predictions 
for the absolute rates for singly-produced pions exceed the 
observed rates in both analyses, however the discrepancy 
is less severe for production by i/, as reported here. 

5. Conclusion 

The single differential cross sections in 7r° momentum 
and angle have been measured for l7r° production by 
charged-current interactions in plastic scintillator (CH). 
The measurements are found to be in agreement with the 
predictions from GENIE, NuWro, and NEUT event gener¬ 
ators. This agreement is interesting because of the approx¬ 
imations needed to make predictions for this channel. The 
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measured cross section in 7r° momentum disfavors the GE¬ 
NIE prediction without FSI effects at low 7r° momentum, 
which is not surprising since FSI is expected for hadrons 
inside nuclei. A decomposition of the FSI effects shows 
that inelastic scattering and CEX reactions are responsi¬ 
ble for the peak at p n ~ 0.2 GeV/c. These contributions 
could be adjusted within external experimental errors to 
achieve even better agreement of calculations with the 7 r° 
momentum data. However, these changes would not be as 
effective for the 7r° polar angle. 

Charged-current single pion production in the few GeV 
region of neutrino energy is an important class of inter¬ 
actions in long-baseline neutrino oscillation experiments. 
This work presents the first detailed measurements of the 
kinematic distributions for single 7r°s produced in charged 
current interactions on carbon. These distributions pro¬ 
vide constraints on antineutrino-induced 7 r° backgrounds 
as will occur in v e appearance experiments. 

Acknowledgments 

This work was supported by the Fermi National Accel¬ 
erator Laboratory under US Department of Energy con¬ 
tract No. DE-AC02-07CH11359 which included the MIN- 
ERvA construction project. Construction support also 
was granted by the United States National Science Foun¬ 
dation under Award PHY-0619727 and by the University 
of Rochester. Support for participating scientists was pro¬ 
vided by NSF and DOE (USA) by CAPES and CNPq 
(Brazil), by CoNaCyT (Mexico), by CONICYT (Chile), 
by CONCYTEC, DGI-PUCP and IDI/IGI-UNI (Peru), by 
Latin American Center for Physics (CLAF) and by RAS 
and the Russian Ministry of Education and Science (Rus¬ 
sia). We thank the MINOS Collaboration for use of its 
near detector data. Finally, we thank the staff of Fermilab 
for support of the beamline and the detector. 

References 

References 

[1] L. Alvarez-Ruso, Y. Hayato, J. Nieves, Progress and open ques¬ 
tions in the physics of neutrino cross sections at intermediate 
energies, New J. Phys. 16 (2014) 075015. arXiv: 1403.2673, 
doi:10.1088/1367-2630/16/7/075015. 

[2] D. S. Ayres, et al., NOvA Proposal to Build a 30 Kiloton Off- 
Axis Detector to Study Neutrino Oscillations in the Fermilab 
NuMI Beamline, arXiv:hep-ex/0503053. 

[3] T. Akiri, et al., The 2010 Interim Report of the Long-Baseline 
Neutrino Experiment Collaboration Physics Working Groups, 
arXiv:1110.6249. 

[4] A. A. Aguilar-Arevalo, et al., Measurement of muon neu¬ 
trino quasi-elastic scattering on carbon, Phys. Rev. Lett. 100 
(2008) 032301. arXiv:0706.0926, doi:10.1103/PhysRevLett. 
100.032301. 

[5] V. Lyubushkin, et al., A Study of quasi-elastic muon neu¬ 
trino and antineutrino scattering in the NOMAD experiment, 
Eur.Phys.J. C63 (2009) 355—381. arXiv:0812.4543, doi: 10. 
1140/epjc/sl0052-009-1113-0. 


[6] L. Fields, J. Chvojka, et al., Measurement of Muon Antineutrino 
Quasi-Elastic Scattering on a Hydrocarbon Target at E u ~ 3.5 
GeV, Phys. Rev. Lett. Ill (2013) 022501. arXiv: 1305.2234. 

[7] A. A. Aguilar-Arevalo, et al., Measurement of neutrino-induced 
charged-current charged pion production cross sections on min¬ 
eral oil at E u ~ 1 GeV, Phys. Rev. D 83 (2011) 052007. 
doi:10.1103/PhysRevD.83.052007. 

[8] O. Lalakulich and U. Mosel, Pion production in the Mini- 
B 00 NE experiment, Phys. Rev. C 87 (2013) 014602. doi: 
10.1103/PhysRevC.87.014602. 

[9] E. Hernandez, J. Nieves, and M.J. Vicente Vacas, Single n pro¬ 
duction in neutrino nucleus scattering, Phys. Rev. D 87 (2013) 
113009. doi:10.1103/PhysRevD.87.113009. 

[10] P. Rodrigues, Comparing pion production models to Mini- 
B 00 NE dataarXiv: 1402.4709. 

[11] J. Y. Yu, E. Paschos, I. Schienbein, Comparison of the ANP 
model with the data for neutrino induced single pion pro¬ 
duction from the MiniBooNE and MINERVA experiments, 
Phys.Rev. D91 (5) (2015) 054038. arXiv: 1411.6637, doi: 
10.1103/PhysRevD.91.054038. 

[12] H. J. Grabosch, et al., Z. Phys. C 41 (1989) 527. 

[13] D. Allasia, et al., Nucl. Phys. B 343 (1990) 285. 

[14] S. Barish, et al., Phys. Rev. D 19 (1979) 2521. 

[15] G. Radecky, et al., Phys. Rev. D 25 (1982) 1161. 

[16] T. Kitagaki, et al., Phys. Rev. D 34 (1986) 2554. 

[17] A. A. Aguilar-Arevalo, et al., The MiniBooNE Detector, Nucl. 
Instrum. Meth. A 599 (2009) 28-46. arXiv:0806.4201, doi: 
10.1016/j.nima.2008.10.028. 

[18] C. Mariani, Neutral pion cross section measurement at K2K, 
AIP Conf.Proc. 967 (2007) 174-178. doi: 10.1063/1.2834471. 

[19] Y. Kurimoto, et al., Measurement of Inclusive Neutral Cur¬ 
rent Neutral 7 r° Production on Carbon in a Few-GeV Neu¬ 
trino Beam, Phys.Rev. D81 (2010) 033004. arXiv:0910.5768, 
doi:10.1103/PhysRevD.81.033004. 

[20] B. Eberly, et al., Charged Pion Production in Interactions 
on Hydrocarbon at (E u )= 4.0 GeVarXiv: 1406.6415. 

[21] A. Aguilar-Arevalo, et al., Phys. Rev. D 83 (2011) 052009. 

[22] D. Rein, L. M. Sehgal, Neutrino Excitation of Baryon Reso¬ 
nances and Single Pion Production, Annals Phys. 133 (1981) 
79-153. doi:10.1016/0003-4916(81)90242-6. 

[23] K. Abe, et al., The T2K Experiment, Nucl. Instrum. Meth. A 
659 (2011) 106-135. arXiv:1106.1238, doi:10.1016/j.nima. 
2011.06.067. 

[24] A. A. Aguilar-Arevalo, et al., Event Excess in the Mini¬ 
BooNE Search for —> v e Oscillations, Phys. Rev. Lett. 105 

(2010) 181801. arXiv:1007.1150, doi:10.1103/PhysRevLett. 
105.181801. 

[25] C. Andreopoulos, et al., The GENIE neutrino Monte Carlo gen¬ 
erator, Nucl. Instrum. Meth. A 614 (1) (2010) 87 - 104, Program 
version 2.6.2 used here, doi : 10.1016/j .nima.2009.12.009. 

[26] C. Juszczak, J. A. Nowak, J. T. Sobczyk, Simulations from 
a new neutrino event generator, Nucl. Phys. Proc. Suppl. 
159 (2006) 211-216. arXiv:hep-ph/0512365, doi:10.1016/j. 
nuclphysbps.2006.08.069. 

[27] T. Golan, J. Sobczyk, J. Zmuda, NuWro: the Wroclaw Monte 
Carlo Generator of Neutrino Interactions, Nucl. Phys. Proc. 
Suppl. 229-232 (2012) 499. doi: 10.1016/j .nuclphysbps .2012. 
09.136. 

[28] Y. Hayato, A neutrino interaction simulation program library 
NEUT, Acta Phys.Polon. B40 (2009) 2477-2489. 

[29] K. Anderson, et al., The NuMI Facility Technical Design Re¬ 
port, Fermilab-design-1998-01 (1998). 

[30] S. Agostinelli, et al., Geant4 simulation toolkit, Nucl. In¬ 
strum. Meth. A 506 (3) (2003) 250 - 303. doi: 10.1016/ 
S0168-9002(03)01368-8. 

[31] J. Allison, et al., Geant4 developments and applications, Nu¬ 
clear Science, IEEE Transactions on 53 (1) (2006) 270—278. 
doi:10.1109/TNS.2006.869826. 

[32] C. Alt, et al., Inclusive production of charged pions in p+C 
collisions at 158-GeV/c beam momentum, Eur. Phys. J. C49 
(2007) 897—917. arXiv:hep-ex/0606028, doi:10.1140/epjc/ 


7 



s10052-006-0165-7. 

[33] D. Barton, G. Brandenburg, W. Busza, T. Dobrowolski, 
J. Friedman, et al., Experimental Study of the a-Dependence of 
Inclusive Hadron Fragmentation, Phys.Rev. D27 (1983) 2580. 
doi:10.1103/PhysRevD.27.2580. 

[34] A. V. Lebedev, Ratio of pion kaon production in proton carbon 
interactions, FERMILAB-THESIS-2007-76. 

[35] L. Aliaga, et al., Design, Calibration, and Performance of the 
MINERvA Detector, Nucl. Instrum. Meth. A 743 (2013) 130. 
arXiv:1305.5199. 

[36] D. G. Michael, et al., The Magnetized steel and scintillator 
calorimeters of the MINOS experiment, Nucl. Instrum. Meth. 
A 596 (2008) 190-228. arXiv:0805.3170, doi:10.1016/j.nima. 
2008.08.003. 

[37] L. Aliaga, et al., MINERvA neutrino detector response mea¬ 
sured with test beam data, Nucl.Instrum.Meth. A789 (2015) 
28-42. arXiv:1501.06431, doi:10.1016/j.nima.2015.04.003. 

[38] J. Beringer, et al., Review of Particle Physics (RPP), Phys. Rev. 
D 86 (2012) 010001. doi:10.1103/PhysRevD.86.010001. 

[39] A. Bodek, I. Park, U.-K. Yang, Improved low Q 2 model for 
neutrino and electron nucleon cross sections in few GeV region, 
Nucl. Phys. Proc. Suppl. 139 (2005) 113-118. arXiv :hep-ph/ 
0411202, doi:10.1016/j.nuclphysbps.2004.11.208. 

[40] G. Radecky, et al., Study of Single Pion Production by Weak 
Charged Currents in Low-energy Neutrino d Interactions, Phys. 
Rev. D 25 (1982) 1161-1173. 

[41] T. Kitagaki, et al., Study of vd —> p~pp s and ud —> 
p~ (1232)n s using the BNL 7-foot deuterium filled bub¬ 
ble chamber, Phys. Rev. D 42 (1990) 1331-1338. doi: 10.1103/ 
PhysRevD.42.1331. 

[42] T. Lee, R. Redwine, Pion nucleus interactions, Ann. Rev. Nucl. 
Part. Sci. 52 (2002) 23-63. 

[43] D. Ashery, et al., True Absorption and Scattering of Pions on 
Nuclei, Phys. Rev. C 23 (1981) 2173-2185. 

[44] C. Anderson, et al., The ArgoNeuT detector in the NuMI low- 
energy beam line at Fermilab, Journal of Instrumentation 7 
(2012) 10019. 

[45] R. J. Barlow, C. Beeston, Fitting using finite Monte Carlo 
samples, Comput.Phys.Commun. 77 (1993) 219—228, as imple¬ 
mented in TFractionFitter, ROOT framework, doi: 10.1016/ 
0010-4655(93)90005-W. 

[46] G. D’Agostini, A multidimensional unfolding method based on 
Bayes’ theorem, Nucl. Instrum. Meth. A 362 (1995) 487-498. 
doi:10.1016/0168-9002(95)00274-X. 

[47] W. P. Abfalterer, et al., Measurement of neutron total cross- 
sections up to 560-MeV, Phys. Rev. C 63 (2001) 044608. doi: 
10.1103/PhysRevC.63.044608. 

[48] W. Schimmerling, et al., Neutron-nucleus total and inelastic 
cross-sections - 900 to 2600 MeV/c, Phys. Rev. C 7 (1973) 248- 
262. doi:10.1103/PhysRevC.7.248. 

[49] R. G. P. Voss, R. Wilson, Proc. Roy. Soc. A236 (1956) 41. 

[50] I. Slypen, V. Corcalciuc, J. P. Meulders, Proton and deuteron 
production in neutron-induced reactions on carbon at E(n) = 
42.5, 62.7, and 72.8 MeV, Phys. Rev. C 51 (1995) 1303-1311. 
doi:10.1103/PhysRevC.51.1303. 

[51] J. Franz, et al., Neutron induced production of protons, 
deuterons and tritons on copper and bismuth, Nucl. Phys. A510 
(1990) 774-802. doi:10.1016/0375-9474(90)90360-X. 

[52] U. Tippawan, et al., Light-ion production in the interaction of 
96-MeV neutrons with carbon, Phys. Rev. C 79 (2009) 064611. 
arXiv:0812.0701, doi:10.1103/PhysRevC.79.064611. 

[53] R. Bevilacqua, et al., Light-ion production from O, Si, Fe and Bi 
induced by 175 MeV quasi-monoenergetic neutronarXiv: 1303. 
4637. 

[54] C. I. Zanelli, et al., Total non-elastic cross sections of neutrons 
on C, O, Ca, and Fe at 40.3 and 50.4 MeV, Phys. Rev. C 23 
(1981) 1015-1022. doi:10.1103/PhysRevC.23.1015. 

[55] T. Bowles, D. Geesaman, R. Holt, H. Jackson, J. Julien, et al., 
Inclusive (tt^, 7r°) Reactions in Nuclei, Phys.Rev. C23 (1981) 
439-447. doi:10.1103/PhysRevC.23.439. 

[56] M. Jones, R. Ransome, V. Cupps, R. Fergerson, C. Mor¬ 


ris, et al., Pion absorption above the Delta (1232) resonance, 
Phys.Rev. C48 (1993) 2800-2817. doi:10.1103/PhysRevC.48. 
2800. 

[57] K. Graczyk, D. Kielczewska, P. Przewlocki, J. Sobczyk, 
axial form factor from bubble chamber experiments, Phys. Rev. 
D 80 (2009) 093001. arXiv:0908.2175, doi:10.1103/PhysRevD. 
80.093001. 

[58] L. Salcedo, E. Oset, M. Vicente-Vacas, C. Garcia-Recio, Com¬ 
puter simulation of inclusive pion nuclear reactions, Nucl. Phys. 
A484 (1988) 557. doi:10.1016/0375-9474(88)90310-7. 



6. Appendix 

This Appendix contains additional tables that are ref¬ 
erenced in the article. 


Table 1: Flux-averaged differential cross section in 7r° momentum, 
dcr/dp^o (10 -40 cm 2 /nucleon/(GeV/c), for l7r° production with sta¬ 
tistical (stat), systematic (sys), and total (tot) uncertainties. 


Pw° (GeV/c) 

da / dp *o 

stat(%) 

sys(%) 

tot(%) 

0.00 - 0.08 

3.75 

41 

33 

53 

0.08 - 0.14 

22.60 

14 

21 

25 

0.14 - 0.20 

27.75 

10 

18 

21 

0.20 - 0.28 

17.92 

11 

21 

24 

0.28 - 0.36 

14.26 

11 

20 

23 

0.36 - 0.45 

12.77 

10 

20 

22 

0.45 - 0.55 

8.65 

11 

20 

23 

0.55 - 0.65 

5.86 

13 

22 

26 

0.65 - 0.75 

4.67 

12 

20 

24 

0.75 - 1.00 

2.25 

11 

18 

21 

1.00 - 1.50 

1.27 

15 

21 

26 


Table 2: Flux-averaged differential cross section in 7r° angle, 

da/dO n o (10 _42 cm 2 /nucleon/deg.), for l7r° production with statis¬ 
tical (stat), systematic (sys), and total (tot) uncertainties. 


9*o (deg.) 

da/d9*o 

stat(%) 

sys(%) 

tot(%) 

0 - 10 

3.97 

34 

19 

40 

10 - 20 

14.25 

15 

15 

22 

20 - 30 

12.88 

13 

15 

21 

30 - 40 

14.06 

12 

15 

19 

40 - 50 

12.70 

12 

15 

19 

50 - 60 

10.32 

14 

15 

21 

60 - 70 

8.14 

17 

15 

23 

70 - 80 

8.95 

17 

16 

23 

80 - 90 

5.11 

21 

18 

28 

90 - 120 

4.73 

16 

18 

24 

120- 180 

1.75 

20 

20 

28 
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SUPPLEMENTAL MATERIAL 


This supplemental material contains additional tables 
that are referenced in the article. 


TABLE I: Flux-averaged differential cross section in 7r° mo¬ 
mentum, da/dp^o (10 _40 cm 2 /nucleon/(GeV/c), for l7r° pro¬ 
duction with statistical (stat), systematic (sys), and total 
(tot) uncertainties. 


TV o (GeV/c) 

da/dp n o 

stat(%) 

sys(%) 

tot(%) 

0.00 - 0.08 

3.75 

41 

33 

53 

0.08 - 0.14 

22.60 

14 

21 

25 

0.14 - 0.20 

27.75 

10 

18 

21 

0.20 - 0.28 

17.92 

11 

21 

24 

0.28 - 0.36 

14.26 

11 

20 

23 

0.36 - 0.45 

12.77 

10 

20 

22 

0.45 - 0.55 

8.65 

11 

20 

23 

0.55 - 0.65 

5.86 

13 

22 

26 

0.65 - 0.75 

4.67 

12 

20 

24 

0.75 - 1.00 

2.25 

11 

18 

21 

1.00 - 1.50 

1.27 

15 

21 

26 


TABLE II: Flux-averaged differential cross section in 7r° 
angle, da/dd^o (lCU 42 cm 2 /nucleon/deg.), for l7r° production 
with statistical (stat), systematic (sys), and total (tot) uncer¬ 
tainties. 


9 w o (deg.) 

da/dd n o 

stat(%) 

sys(%) 

tot(%) 

0 - 10 

3.97 

34 

19 

40 

10 - 20 

14.25 

15 

15 

22 

20 - 30 

12.88 

13 

15 

21 

30 - 40 

14.06 

12 

15 

19 

40 - 50 

12.70 

12 

15 

19 

50 - 60 

10.32 

14 

15 

21 

60 - 70 

8.14 

17 

15 

23 

70 - 80 

8.95 

17 

16 

23 

80 - 90 

5.11 

21 

18 

28 

90 - 120 

4.73 

16 

18 

24 

120- 180 

1.75 

20 

20 

28 



